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ABSTRACT 

Values  of  four  parameterization  constants  ror  che  duik, 
second-order  closure  turbulence  model  of  Garwooa  (1976)  are 
determined  using  ocean  data  from  the  Nortn  Paciric  experi- 
ment (NORPAX)  Trans-Pacific  ship  of  opportunity  program 
(TRANSPAC)  ana  atmospneric  forcing  from  the  Fleet  Numerical 
Oceanography  Center  (FNOC)  mooei.  Tne  annual  cycle  or  sea 
surface  temperature  is  simulatea  ror  eacn  of  tne  years 
1976-1978.  The  difference  oetween  tne  simulatea  ana  tne  an- 
alyzed monthly  change  of  sea  surface  temperature  is  mini- 
mized to  determine  the  optimal  values  of  trie  moaei  con- 
stants. The  moael  is  relatively  insensitive  to  cnanges  in 
the  constants  for  this  set  of  aata.  However  simulations  us- 
ing the  values  aeterminea  in  tms  stuay  yieia  results  wnicn 
are  quite  different  from  those  obtainea  using  vaiues  ror  tne 
model  constants  aeterminea  from  earlier  stuaies  witn  airrer- 
ent  initialization  ana  forcing  aata. 

Three  statistical  techniques  are  usea  to  determine  tne 
sensitivity  of  the  moael  to  cnanges  of  tne  moaei  constants. 
Two  techniques,  bootstrapping  ana  pooiea  permutation  proce- 
dure, yield  similar  conclusions.  The  tnira,  a  stanaara  stu- 
dent's t-test,  provea  to  oe  overiy  sensitive.  This  leaas  co 
the  conclusion  that  Gaussian  statistics  are  not  appropriate 
to  studies,  such  as  this  one,  for  wnicn  tne  aegrees  of  rree- 
dom  are  not  well  defined. 

The  formulation  ror  tne  aosorption  or  solar  raaiation 
witn  depth,  in  tne  Garwooa  moael,  is  a  simplification  of  a 
two  exponent  model.  It  aaequateiy  reproauces  tne  character- 
istics of  more  complex  absorption  moaeis.  Optimization  of 
the  constants  of  the  absorption  model  improves  tne  simula- 
tions of  the  summer  mixea  layer  ana  the  annual  cycle  or  cne 
mixea  layer. 
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1.   INTRODUCTION 

This  paper  aaaresses  three  interrelated  topics  wnicn 
are  fundamental  to  simulating  the  thermodynamics  or  tne  up- 
per ocean  on  aiurnai  to  seasonal  time  scales.  Tne  rirst 
topic  is  the  determination  tne  mooel  constants  rcr  tne  sec- 
ond-order closure  bulk  mooel  of  turouient  oceanic  oounoary 
layers  aevelopeo  by  Garwooa  (1976,  1977).  Tne  constants  are 
determineo  for  forcing  aeriveo  from  the  Fleet  Numerical  Gce- 
anographic  Center  (FNQC)  atmospheric  primitive  equation  moo- 
el and  ocean  thermal  structure  information  obtained  from  the 
Nortn  Pacific  Experiment  (NGRPAX)  Trans-Pacific  snip  of  op- 
portunity XBT  ( TRANS PAC)  analysis.  Seconuiy,  tnree  statis- 
tical tests  are  usea  to  examine  tne  variations  in  mooel 
skill  as  a  runction  of  tne  parameterization  constants  ano 
the  solar  raaiation  absorption  mooel.  Finally,  tne  model 
for  the  absorption  of  solar  radiation  in  tne  upper  ocean, 
which  is  used  in  tne  Garwooa  moael,  is  examined  ano  compared 
with  several  otner  moaels. 

Turbulence  closure  is  necessitatea  uy  tne  process  or 
Reynold's  averaging  whicn  generates  correlations  (moments) 
of  tne  fluctuating  variaoies,  resulting  in  more  unknowns 
than  equations.  The  proDlem  can  not  De  resoiveu  at  a  funda- 
mental level  with  adoitional  equations  since  each  new  equa- 
tion for  a  correlation  mtroouces  aoaitionai  nigner  oroer 
moments.  The  system  of  equations  must  oe  cioseo  oy  speciry- 
ing  the  hignest  oroer  of  the  terms  to  oe  retaineu  ano  param- 
eterizing these  correlations  in  terms  of  iower  oroer  moments 
ano  mean  values.  One  may  tnen  characterize  the  various  tur- 
bulence mooeis  accoromg  to  tne  oroer  at  wnicn  this  closure 
is  mace  ano  the  particular  parameterizations  usea.  Tne  Gar- 
wooa (1976,1977)  moael  aiscussea  in  the  paper  employs  a  sec- 
ona-order  closure  scneme.  Thus,  tne  secona-oroer  moments 
(triple  correlations)  are  parameterizea  in  terms  of  first- 
order  moments  (autocorrelations  and  cross  correlations)  ana 
zero-order  moments  (mean  values) . 

The  parameterization  or  tne  nignest  order  momencs  gen- 
erates constants  of  proportionality   oetween  tne  nignest  or- 
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der  moments  and  functions  of  tne  lower  oraer  moments.  These 
constants  must  be  determined  experimentally  berore  the  mouex 
can  be  used  to  simulate  accurately  the  fluid  flows  or  inter- 
est. The  goal  in  turbulence  modeling  is  to  formulate  tne 
parameterizations  such  that  tney  are  "universal"  or  applica- 
ble to  a  wide  range  of  geophysical  flows.  Unfortunately,  no 
parameterization  can  be  completely  universal.  Parameteriza- 
tions derivea  for  bounaary  layer  flows,  for  example,  can  not 
be  expected  to  work  well  for  free  mixing  layers,  jets  or 
wakes,  since  the  physics  ana  tne  oounaary  conditions  are 
quite  different.  Further,  for  a  given  parameterization,  the 
optimal  values  of  the  constants  may  also  oe  expecteu  to  an- 
fer  for  different  oceanic  ana  atmospheric  regimes  since  all 
possible  pnysical  processes  are  not  inciuaea  in  the  moael. 
In  adaition,  the  space  or  time  resolution  of  tne  aata  usee 
may  not  be  adequate.  The  constants  for  this  moael  nave  pre- 
viously been  estimatea  form  laboratory  aata  (Garwooa,  1976) 
and  ocean  weather  ship  observations  at  Ocean  Station  Papa 
(50%,  14oW)  (Garwood  ana  Aaamec,  1982)  .  Although  tnese  ear- 
lier values  were  not  optimal  aeterminations,  tney  will  oe 
compared  with  the  reults  ootainea  nere  using  FNOC  ana 
TRANS PAC  aata. 

The  absorption  of  solar  raaiation  in  the  upper  ocean 
significantly  inriuences  tne  aiurnai  cycle  of  the  mixeo 
layer,  the  summer  minimum  mixea  layer  deptn  ana  associatea 
temperature  maximum  ana  the  spring  transition.  The  physics 
of  the  absorption  of  solar  raaiation  is  v/ell  known;  nowever, 
the  distribution  of  yellow  matter  ana  particulate  matter, 
which  significantly  affects  tne  process,  is  not.  Further, 
an  elaborate  banawidth  aepenaent  moael  of  solar  raaiation 
absorption  is  very  computationally  expensive.  Thus  it  is 
necessary  to  develop  an  accurate  out  simpiifiea  parameteri- 
zation of  the  absorption  process.  The  formulation  usea  in 
this  moael  will  be  comparea  with  other  methoas,  inducing  an 
empirical  moael  (Zanevela  ana  Spinraa,  1980)  wnich  oest  rep- 
resents the  limited  aata  sets  that  are  avaiiaoie. 
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The  choice  ana  application  of  statistics  wmcn  are  ap- 
propriate to  this  type  of  stuay  is  a  non-triviai  aspect  or 
the  work.  Three  relatively  simple  statistics  nave  oeen  cho- 
sen: one  Gaussian  statistic  ana  two  non-Gaussian  statistics 
which  are  capable  of  generating  their  own  prooaDixity  at- 
tribution functions  (paf's).  It  is  nopea  that  tnis  latter 
approach  will  reduce  the  proolems  causea  Dy  using  sample 
populations  which  do  not  necessarily  nave  Gaussian  paf's  ana 
whose  members  are  not  inaepenuent. 

The  Garwooa  moael  is  reviewea  m  section  2.  Tne  ocean 
thermal  structure,  surface  forcing  ana  tne  aetails  of  tne 
numerical  simulations  for  this  stuay  are  aiscussea  in  sec- 
tion 3.  The  results  of  the  numerical  experiments  are  aa- 
dressed  in  section  4  ana  a  summary  is  proviaea  in  section  5. 


2.   DESCRIPTION  OF  MODEL 

Since  the  moael  has  been  descrioea  in  aetail  Dy  Garwooa 
(1976  ana  1977) ,  we  will  only  briefly  review  the  equations 
to  indicate  the  moael  constants  which  must  oe  aeterminea. 

Garwooa  integrated  the  equations  of  motion  ana  trie 
buoyancy  equation  through  the  mixea  layer  to  yield 

3<U> 


n-sr -  f  h<V>  =  -uw(0)  +  uw(-n)  , 

dt 
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The  vertical  mean  value  of  any  vanaole,  X,  is  uerineu  as 
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=  -   X(z)az  , 


(2.3) 


where  h  is  the  mixed  layer  aeptn  ana  z  is  zero  at  tne  sea 
surface  and  positive  upward.  The  oounaary  conaitions  at  z  = 
0  ana  z  =  -h  are 
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Tx  and  11  are  the  x  ana  y  components  of  the  surface  stress. 
Qv  is  the  net  surface  heat  flux  (total  heat  flux,  QT,  less 
solar  radiation,  Qft) .  S  is  the  surface  salinity  ana  (E-P) 
is  the  evaporation  less  precipitation  rate.  _/\_,  the  Heav- 
iside  step  function,  is  aefinea  as 


a  dh 

_A_  =  0  for  —  <  0 
dt 


ana    _/\_  =  1  ror  —  >  0  . 

at 


The  changes  in  mean  velocity  components  ana  ouoyancy  tnrougn 
the  entrainment  zone  are 

Au  =  U(-h)  -  U(-h  -  S)    , 

Av  =  V(-h)  -  V(-n  -  <J)  , 

Ab  =  B(-h)  -  B(-n  -  S)    . 
/\u ,      /\V  and  Ab   are  the  cnanges   in  tne  mean  properties 
across  the  entrainment   zone  ana  ah/at  is  tne   rate  of  mixea 
layer  deepening. 

The  turbulence   equations  integratea  througn   the  mixea 
layer  are 
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m, 


Integrating  the  total  turbulent  kinetic  energy  (TKE)  tnrougn 
the  entrainment  zone  proviaes  the  final  equation   neeaea  to 
close  the  system  of  equations, 
2m¥(TAA<e~> 


dh 

dt   hAB  +  2<e> 
<e>  is  the  vertically  averagea  TKE 


(2.6) 
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—    <u*>  +  <v*>  +  <w3»> 
2~     2  ~ 

The  left  hand  siae  of  (2.4)   is   the  local  rate  or  change  of 
the  vertically  averaged  horizontal  TKE  in  tne  raixea  layer. 

The  time  scale  for  the  terms  on  the  left  nanu  siue  of 
(2.3)  and  (2.4)  is  approximately  100  seconas  for  tne  oceanic 
planetary  ooundary  layer  under  normal  circumstances  (ae 
Szoeke  and  Rhines,  1976).  Since  we  are  interestea  in  montn- 
ly  ana  longer  time  scales  in  tnis  stuay,  tnese  terms  will  oe 
neglected.  The  first  term  on  tne  rignt  nana  siae  of  (2.3) 
is  a  parameterization  of  the  surface  wina  stress  actmy  on 
the  shear  of  the  mean  current,  tne  proauction  aue  to  creat- 
ing surface  waves  ana  the  transport  or  norizontai  TKE.  u^ 
is  the  ocean  friction  velocity,  given  oy  u*  =  iT|/o  ana 
used  as  a  scale  for  the  rms  velocity  fluctuations.  Trie  sec- 
ond term  on  the  right  hand  siae  of  (2.3)  is  tne  momentum 
flux  associatea  with  the  entrainment  of  nonturouient  fluid 
at  the  base  of  the  mixea  layer.  The  momentum  rlux  aue  to 
internal  waves  is  neglectea.  Tne  thira  term  is  a  parameter- 
ization of  the  pressure  rate-of-strain  terms  wnicn  act  to 
reduce  the  anisotropy  of  tne  turouient  flow  (Rotta,1951; 
Lumley  and  Khajeh-Nouri ,1974) .  The  last  term  is  a  parame- 
terization of  the  aissipation  rate.  The  aissipation  rate 
can  be  estimatea  from  the  rate  at  wnicn  tne  large  scale  ea- 
dies  supply  energy  to  the  smaller  scales.  As  aiscussea  in 
Tennekes  ana  Lumley  (1972)  ,  the  dissipation  rate  is 

OU:  du; 

E 


C       flu-  du;      <e> 

=  1  r  T-^-r-  az  =  —  ' 

JK  dxjAxj       Tfc 


where  T  is  the  dissipation  time  scale.  There  are  two  time 
scales  for  the  oceanic  mixed  layer:  tne  time  scale  imposed 
by  planetary  rotation,  T^^l/f  ana  the  time  scale  formed 
from  the  turbulent  velocity  ana  the  length  scale  of  tne 
large  scale  turouient  flow,  T-***  h/u^..  The  first  scale  is 
the  characteristic  time  scale  of  vortex  stretcning  aue  co 
the  planetary  rotation.  The  secona  is  tne  time  re^uirea  ror 
a  typical  large  scale  eaay  to  overturn.  A  reature  or  Gar- 
wooa's  model  is  tnat  Doth  of   tnese  time  scales  are  mcorpo- 
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rated   into  the  parameterization  of   the  aissipation   tii&e 

scale  by  specifying  that  T£  is   scaled  with  the  smaller  or  T^ 

or  T_.   To  incorporate  both  time  scales  the  aissipation  time 

scale  is  defined  as 
-»    -i    -I 

T£  =  TR  +  TT  . 

In  (2.5)  ,  the  first  two  terms  on  the  rignt  nana  siae 
are  the  surface  and  entrainment  buoyancy  fluxes.  Tne  tnira 
term  is  the  absorption  of  solar  raaiation  witn  aeptn.  Tne 
remaining  terms  are  identical  to  those  in  (2.4).  (2.1) 
through  (2.6)  form  a  complete  set  or  equations  for  trie  de- 
pendent variables  <U>,  <V> ,  <B>,  <u3»  +  v*»>,  <w^>  ana  n. 

The  depth-depenaent  absorption  of  solar  raaiation  is 
particularly  important  in  shallow,  summer  mixea  layers.  In 
these  cases  the  amount  of  solar  raaiation  aosoroea  in  tne 
mixea  layer  significantly  affects  tne  mixea  layer  tempera- 
ture, while  the  absorption  oeiow  the  mixea  layer  arrects  tne 
entrainment  rate  and,  tnerefore,  the  mixea  layer  aeptn.  Tne 
flux  of  solar  radiation  in  the  ocean  ean  oe  written  as 

where  h   is  the  wavelength  ana  W(^)  is  tne  optical  water  type 
(Jerlov,   1976).   Various  parameterizations   (see  TaDie  1) 
have  been  proposea  for  the  aosorption  of  solar  raaiation. 

TABLE  1 
Parameterizations  of  solar  radiation  in  tne  ocean 

Parameterizations  Source 

QJz)    =  Qpe**  Denman,    1973 

Qs(z)    =   ^)(r,  eX*z    +    raey**)  Paulson   ana   Simpson,    1977 

Q^z)    =   QU)Z  r;  eYii  Konao,    197  9 

Qs(z)    =   <^eyi(l   -   a   tan"' (ox))  Zaneveia   ana  Spmraa,    19b0 

Q^z)    =  QfaHl-r)  |z|    <  1  Garwooa,    1976 

=   Qtfre*2  |z|    >    1 

Q$ut)is   the  solar  radiation  at  the  sea  surface,  tne  0 '  s  are  the 
extinction  coefficients  ana  the  r's  are  constants  (0  <  r  <  1). 

The  first   three  of   these  parameterizations   assume  that  a 
significant  fraction   of  the  solar  raaiation   aosorption  can 
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be  modeled  by  using  one,  two  or  many  wavelength  Danes,  re- 
spectively. The  fourth  parameterization  was  ootainea  empir- 
ically using  an  extensive  data  set  of  absorptance  ana  trans- 
mittance  measurements  (Spinrad  et.  ai. ,  1979).  Tne  con- 
stants in  the  fourth  formula  in  Taole  1  aepena  on  tne  water 
type.  For  this  study  the  appropriate  water  type  is  Type  II 
(Zaneveld,  personal  communication),  ana  the  constants  are  X 
=  0.067  m/s,  a  =  0.4158,  ana  b  =  3.9865  m/s.  Simpson  ana 
Dickey  (1981)  comparea  the  first  four  parameterizations  in 
Table  1  ana  founa  that  using  a  single  exponential  mouei  is 
acceptable  only  at  very  high  wina  speeas  (U,  >  20  m/s). 
They  also  founa  that  the  empirical  formulation  of  Zaneveia 
and  Spinrad  (1980)  and  the  aouble  exponential  moaei  or  Simp- 
son ana  Paulson  (1977)  proauce  nearly  laentical  temperature 
profiles,  and  that  the  complexity  of  tne  spectral  aecomposi- 
tion  form  is  not  necessary  in  mixea  layer  moaeiing. 

The  fifth  parameterization,  usea  in  tne  Garwooa  moael, 
is  a  douole  exponential  moael  for  wnicn  one  of  tne  extinc- 
tion depths  approaches  zero.  A  fixea  fraction  (1-r)  of  Q 
is  absorbed  in  the  first  meter  ana  the  remaining  flux  is  aD- 
sorbed  exponentially  for  |z|  >  lm.  Values  of  r  ana  o  must 
be  specified.  This  parameterization  will  oe  comparea  witn 
Zanevela's  empirical  formulation. 

In  (2.4)    to  (2.6),   m,   througn  m-  are   tne  constants 
which  must  be   determined  from   iaooratory  anu  geophysical 
data.   Of  these  only  four  are  inaepenaent,  ana  only  tnree  (m- 
and  p3  =  ms/mf  ) 

are  critical  to  moaeiing  tne  annual  cycle  of  the  upper 
ocean  mixed  layer.  In  tnis  stuay,  we  use  neat  ana  momentum 
fluxes  derived  from  the  FNOC  atraospneric  fields  ana  me 
TRANSPAC  ocean  data  over  a  range  of  latitudes  ana  longituaes 
to  estimate  the  optimal  values  for  m _ ,  p  =  m^/n^  ,  r  ana  o. 

3.   DESCRIPTION  OF  NUMERICAL  EXPERIMENTS 
Choosing  an  annual  time  periou   for  this  stuay  requires 
that  the  moael  be  able  to  simulate  four  signiricant  pnenome- 
na  each  of  which  requires   aifferent  responses  rrom  tne  ruoa- 


-  7  - 


el.  The  fall  deepening  is  particularly  sensitive  to  m-  (eq. 
2.4),  the  parameterization  of  surface  wind  mixing.  The  win- 
ter minimum  mixed  layer  temperature  is  sensitive  to  p  (eqs. 
2.4,2.5),  the  dissipation  parameterization.  Trie  spring 
transition  aepends  on  the  interaction  of  m  ano  tne  surrace 
heat  fluxes.  Finally,  the  summer  maximum  mixeo  layer  temp- 
erature depenas  on  r  ana  X  in  the  parameterization  of  solar 
raaiation  absorption. 

A  grid  covering  a  large  area  of  the  North  Paciric  rrom 
34  M  to  46  H  in  4  increments  and  from  170  W  to  140  W  in  10 
increments  was  cnosen  for  this  stuuy.  Tnis  area  is  large 
enough  to  encompass  various  oceanic  tnermai  structure  ana 
atmospheric  forcing  regimes  but  aoes  not  inciuae  areas  wmcn 
are  strongly  influencea  oy  eastern  or  western  bounaary  cur- 
rents where  norizontai  processes  (e.g.  aavection  ana  aiver- 
gence)  are  likely  to  be  important.  The  north-south  extent 
of  the  grid  incorporates  only  tne  area  whicn  has  consistent 
data  coverage  throughout  the  entire  year.  Regions  to  the 
north  (south)  of  these  latitudes  have  fewer  aata  points  aur- 
ing  the  winter  (summer)  aue  to  the  latitudinal  shift  of  ship 
tracks  with  the  seasons. 

The  specific  quantity  chosen  to  compare  the  moaei  re- 
suits  with  the  data  is  the  montniy  change  in  sea  surface 
temperature  (SST)  averaged  over  an  annual  time  perioa  ana 
over  the  grid.  The  monthly  change  in  SST  is  more  difficult 
to  simulate  than  the  SST  since  tne  large  scale  nortn-soutn 
SST  gradient  in  the  nortn  Pacific  contains  more  or  tne  vari- 
ance than  that  associatea  with  the  month- to-month  cnange  in 
the  SST  field.    £*:«  4«&nticy  to  ot  minimized  can  ue   written 

aS  MA/ 

E  =  rl  r X £At« - At*J  (3-1) 

where  M  is  the  numoer  of  gria  points,  N  is  the  numDer  of 
months  used  in  the  average,  ana  A.T^  ana  /A^T^are  zne  montniy 
changes  in  the  simulatea  ana  analyzea  SST  respectively. 
This  quantity  was  choosen  rather  tnan  an  cms  difference 
since  it  allows  the  study  of  some  aaaitional  statisitical 
techniques  which  will  oe  reportea  elsewhere. 
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The  ocean  thermal  structure  fields  wnich  are  useu  in 
the  initialization  of  the  model  ana  in  tne  computation  or 
the  model  error  were  ootainea  from  the  NORPAX  TRANSPAC  anal- 
ysis (White  ana  Bernstein,  1979)  .  Tne  optimal  interpolation 
analysis  uses  the  temperature  ooservations  during  a  month  to 
estimate  the  temperatures  at  each  grid  point.  Horizontal 
analyses  are  preparea,  indepenaently ,  at  0 ,  20,  40,  60,  80, 
100,  120,  150,  200,  250,  300,  ana  400  meters.  These  neias 
are  interpolatea  to  5m  aepth  intervals  using  a  comoination 
of  linear  ana  exponential  interpolation  (Eisoerry,  et.  ai. , 
1982) .  The  moael  is  initialized  on  January  15th  ana  inte- 
grated through  Decemoer  15th  for  1976,  1977  ana  1978.  The 
model  generatea  monthly  SST  used  in  (3.1)  is  a  five  aay  av- 
erage centered  on  the  15th  of  each  month. 

The  surface  oounaary  conditions  requirea  ror  the  Gar- 
wooa  moael  are  the  total  heat  flux  (QT) ,  the  solar  raaiation 
(Q  )  and  tne  wina  speed  (W) .  These  variaDles  are  aerivea 
from  fields  generated  oy  the  FNOC  atmospheric  preuiction  ana 
analysis  system.  Instantaneous  values  of  Q  ana  Q  are  ar- 
chivea  at  12-h  intervals  and  W  at  6-n  intervals  on  a  63  x  63 
polar  stereographic  gria  of  the  northern  hemisphere.  Bili- 
near  interpolation  is  usea  to  map  these  neias  onto  a  4  la- 
titude by  10  longituae  gria  cnosen  for  this  stuay.  Tne 
values  are  interpolated  to  1-n  intervals  (Gaiiacner,  1979) 
and  QT  is  corrected  to  remove  a  systematic  long  term  oias 
(Elsberry  et.  al. ,  1982). 

The  proceaure  usea  to  find  tne  oest  fit  of  tne  moael 
simulations  to  the  analyzea  fields  was  to  first  set  r  =  0.5 
ana  o  =  0.1  m"  and  minimize  (3.1)  in  (m_,  p_ )  space  ror 
each  year.  Next,  r  was  nela  constant  ana  X,  m_ ,  ana  p  were 
varied  to  achieve  the  oest  fit  to  the  1977  annual  cycle. 
This  value  of  o  was  comparea  to  tne  optimal  o  f ouna  maepen- 
dently  Oy  comparison  with  Zanevela's  empirical  formula.  In 
tne  comparison  with  Zanevela's  formula  both  r  ana  0  are  var- 
ied. 
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4.   DISCUSSION  OF  RESULTS 

As  a  preliminary  test,   the  1977  annual  cycle  was  usee 

to  determine  the  distrioution  of  the  mociei  error  in  (m— ,   &) 

space.    A  broad  but  physically  reasonaole  range  or  m-   ana  y 

was  tested  using  widely  spacea  values  or  m_ ,  p  .   These  Sim- 
's      i 

ulations  inaicate  a  local  minimum  of  E  in  the  neignoornooa 
of  m3  =  2.0,  p  =  4.0.  The  minima  for  all  three  years  were 
located  by  testing  the  region  centered  on  p =  4.0,  m  =  2.0 
(Figs,  la,  b  ana  c) .  For  each  year  tne  minimum  lies  on  tne 
line  p  =  4.0,  but  the  value  of  m  varies  oetween  1.5  for 
1978  ana  4.6  for  1976. 

4.1   STATISTICS  0_E  VARIATIONS  m   Tflfi  MODEL  CONSTANTS 

The  student's  t-test  is  usea  to  aetermine  wnetner  tne 
variations  in  E  are  statistically  significant.  Tne  null  hy- 
pothesis, He,  is  E(l)  -  E(2)  =  0  wnere  1  ana  2  represent 
any  two  values  of  the  spatial  ana  temporal  average  iaoael  - 
analysis  differences  aefinea  oy  (3.1).  The  nypothesis  is 
accepted  if 

-t  <  t  <  t 

0(/2,m+n-2        o(/2,m+n-2 

where  t  is  the  stuaent  variate,   oC is  tne  significance  level 

and  m  ana  n  are  the  nuniDer  of  points  in  the  samples  (m  =  16, 

n  =10  ;  MxN  =  160,  for  these  stuaies) .    Tne  stuaent  variate 

is  aefined  as 

E(l)  -  E(2) 
t  =  — 

sp  (1/m  +  1/n) 
where  sp,  the  pooiea  variance,  is 

(m-1)  si   +  (n-1)  s2 

sp  = . 

(m+n-2) 

si  and  s2  are  the  variances  of  E(l)  ana  E(2)  respectively. 

For  the  values  in  1976  (Fig.   la),  tne  stuaent's  t-test 

indicates   that  all   non-adjacent  points  are   signir icantly 

different  at   the  90%  confiaence   level  (wnich  will   oe  usea 

for  all  cases) .    Furthermore,  tne  only  aajacent  pairs  of  (ra- 

t    P_)  points  which  ao  not  aiffer  significantly  are  (3.0,  4.0 

and  3.5,  4.0),  (3.2,   5.0  and  3.5,   4.0)   ana  (3.2,   5.0  ana 
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3.2,   4.0).    In  the  1977  simulations,   tne  results  are  mucn 
less  sensitive  to  variations  of  tne  parameters.   Tne  point  m_ 
=  1.5  was   significantly  airferent  trom  the  points   m_  =  2.2 


ana  m_  =  2.4  with  p  =  4.0.  None  of  tne  other  points  were 
significantly  different  for  1977.  In  1978,  with  p  =  4.0, 
the  point  rrij  =  1.8  differs  signir icantiy  from  q_  =  1.5,  1.6, 
ana  3.0  ana  m^  =  2.6  is  significantly  different  from  m_  = 
1.4  ana  m  =  1.5.  If  the  constants  are  set  at  p  =  4.0,  m, 
=  2.2,  the  moael-analysis  aifferences  are  signiricant,  dc 
the  90%  level,  among  the  tnree  years.  Furtnermore,  tne  min- 
imum in  the  parameter  space  for  eacn  year  is  signif icantiy 
different  from  the  minimum  in  eacn  of  tne  otner  two  years. 

The  stuaent's  t-test  assumes  that  tne  samples  are  arawn 
ranaomly  from  a  normal  population.   Thus  the  aata  points  are 
assumed  to  be   inaepenaent  of  eacn  other.     This  assumption 
can  be   tested  by  pooling  the   /^T's  from  two  sets   of  moaei 
simulations.    A  probability  aistrioution  function  (pat)  can 
be  generated  by  ranaomly  arawmg   samples  from  tne  new  popu- 
lation.  This  pdf  can  then  be  usea  to  aetermine  tne  signifi- 
cance levels.    Using  this  methoa  aoes  not  require  that  tne 
samples  have  a  Gaussian  pdf  or  that  tne   points  De  inaepen- 
aent.   An  example  of  this  proceaure,   tne  pool  ana  permuta- 
tion procedure  (PPP)   is  aescrioea  oy  Preisenaorf er  ana  iioo- 
ley  (1982).    First,   tne  two  data   sees  are  pooiea  into  one 
set  where  the  first  m  elements  are  the  first  sample  of  /\t' s 
and  the   following  n  elements   are  tne  secono   sample.    Tne 
next  step  is  to  permute  the  m  +  n  aata  points  to  construct  a 
new  pooled   aata  set.     The  first  m   elements  oecome   a  new 
first  sample  ana  the  following   n  elements  oecome  trie  second 
sample.    The   test  variable   t  may   then  oe   calculatea  ror 
these  two  new  aata  sets.   These  steps  are  now  repeatea  for  a 
total  of   100  random  experiments   yielamg  100  values   or  t. 
The  t  are  then  oraerea  in   increasing  value  witn  t(l)   oeing 
the  smallest  ana   t(100)   tne  largest.    Then   the  values  of 
t(5)  ana  t(95)  now  yield  tne  90%  significance  levels  airect- 
ly,   which   eliminates  the  assumptions   required  to   use  tne 
stuaent's  t-test. 
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A  similar  method  called  'bootstrapping'  can  oe  useG  co 
generate  the  pcif  (Preisenaorf  er ,  personal  communication)  . 
In  bootstrapping,  the  selections  are  maoe  randomly  from  the 
entire  data  set,  whereas  only  the  permutations  of  the  pooiea 
population  are  used  in  the  PPP.  Therefore,  since  repeti- 
tions of  the  data  points  may  enter  into  the  pooieo  sec  using 
bootstrapping,  slightly  broader  reference  distrioutions  may 
be  generated.  Thus,  tne  null  nypothesis  would  teno  to  oe 
accepted  more  often  using  bootstrapping  than  using  tne  PPP. 

The  critical  points  for  accepting  the  nulx  hypothesis 
may  be  quite  different  from  tne  student's  t-test  if  the 
original  data  sets  do  not  have  approximately  Gaussian  pafs. 
If  this  is  the  case,  fewer  runs  will  oe  statistically  dif- 
ferent using  the  oootstrapping  metnoo  and  PPP.  Results  are 
taoulated  in  Taoles  2,  3  and  4  for  cases  in  wmcn  tne  stu- 
dent's t-test  rejected  the  null  nypothesis  at  tne  90%  conri- 
dence  level.  Variations  in  m_  and  p  within  a  given  year 
(Table  2)  do  not  leao  to  any  significant  differences  except 
using  the  PPP  Detween  m.  =  2.0  ana  m_  =  3.5  for  p  =  4.0  in 
1976. 

TABLE  2 
Comparison  of  PPP  ana  Bootstrapping 


Year 

1976 

1976 
1976 
1977 
1978 

Student  variate  ana  90%  confiaence  limits  as  aeterminea  oy  cue 
bootstrapping  metnoa  ana  Oy  PPP  for  the  cases  for  wnicn  tne  stuae 
t-test  rejectea  the  null  nypothesis  at  tne  90%  level  witnm  eacn 
particular  year.  The  entry  marKea  with  an  asteriSK  was  cne  oniy 
case  tor  which  either  test  aeterminea  tnat  the  points  were 
significantly  different. 

Holding  the  parameters  constant  (p  =4.0,  m,  =  2.2)  and 
testing  between  years  (Table  3),  the  PPP  test:  leads  to  tne 
conclusion  that  the  results  for  1976  ano  1977  are  signifi- 
cantly different,  but  tne  oootstrapping  tecnnique  aia  not. 
Both  proceaures  founa  1976  ana  1978  to  oe  significantly  air- 
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(ms,p3) 

( m3  ,  p3  ) 

3.5,4.0) 

3.2,4.0 

3.2,4.0 

1.5,4.0 

1.5,4.0) 

(2.0,4.0 
3.2,5.0 
2.2,4.0 
2.4,4.0 

(3.0,4.0 

t 

PPP 
t(5)     t(95) 

Bootstrapping 
t(5)    tt95) 

-11.49 

-  3.09 

-  8.11 
1.36 

-  5.78 

-11.05      6,66* 
-11.96    10.01 
-12.31      8.02 

-  6.91       8.97 

-  9.46      7.73 

-12. d8      9.72 
-16.27    10.69 
-12.15      9.6  2 

-  8.D8      7 .93 

-  7  .6d      7 .73 

ferent  for  m_  =  3.0,  p  =  4.0  out  not  for  m-  =  2.2,  p_  =  4.0 
(Table  3).  Results  for  1977  anu  1978  are  not  statiscicaixy 
different.  Testing  for  differences  Detween  minima  or  eaon 
year  (Taole  4) ,  both  techniques  founo  that  the  1976  resuxcs 
are  significantly  aifferenc  from  1977  ana  1978,  out  1977  ana 
1978  are  not  significantly  aifferenc 

TABLE  3 
Comparison  of  PPP  ana  Bootstrapping  Between  Years 


Years 

1976  vs  1977 

1976  vs  1978 

1977  vs  1978 
1976  vs  1978 


(rr.3,p3) 

2.2,4.0 
2.2,4.0 
2.2,4.0 
3.0,4.0 


PPP 
t(5)  t(95) 


Bootstrapping 
t(5)  t(35) 


-  8.01 

-  4.88 
-1.83 
-13.16 


-  7.13   9.38*  -  8.82  9.21 

-  7.84  7.86  -  7.75  7.43 
-10.51  7.58  -  7.61  8.81 
-12.08  10.34*  -  9.53  8.93* 


Student  variate  and  90%  confidence  limits  as  determined  by  tne 
bootstrapping  methoa  ana  by  PPP  for  the  cases  ror  wnicn  tne  stuaent's 
t-test  rejectea  the  null  nypotnesis  at  tne  90%  level.  Tne  entries 
marked  with  asterisks  were  the  only  cases  for  wnicn  eitner  test 
determined  that  the  points  were  significantly  aifrerent. 


TABLE  4 
Comparison  of  PPP  ana  Bootstrapping  Between  Minima 


Years 

1976  vs  1977  -17.49 

1976  vs  1978  -  9.55 

1977  vs  1978  2.42 


PPP 
t(5)  t(95) 

6.59  7.47* 
7.42  7.44* 
7.88   6.00 


Bootstrapping 
t(5)   t?95) 

-12.18   5.86* 

-  9.22   6.90* 

-  9.35   6.47 


Student  variate  ana  90%  confiaence  limits  as  aeterminea  by  the 
bootstrapping  methoa  and  oy  PPP  for  the  minima  oetween  eacn  or 
the  three  years.  These  (mf/p  )  minima  are:  1976  (4.6,4.0); 
1977  (2.4,4.0);  ana  1978  1 1  ?  5  ,4  .0 )  . 

The  entries  markea  with  asterisks  were  cases  for  wnicn  the  test 
determined  that  the  points  were  significantly  aifrerent. 


For  most  cases,  the  oootstrapping  metnoa  yieias  smaller 
values  of  t(5%)  ana  larger  values  of  t(95%)  than  aoes  tne 
PPP.  In  aaaition,  both  techniques  yiela  smaller  values  of 
t(5%)  and  larger  values  of  t(95%)  than  the  stuaent's  t-test. 
Again,  this  implies  that  tne  oootstrapping  metnoa  ana  tne 
PPP  will  fina  two  cases  to  be  significantly  aifrerent  less 
often  than  will  the  stuaent's  t-test.  This  large  airrerence 
in  tne  critical  values  usea  for  the  aifrerent  types  of  tests 
inaicates  that  the  assumption  of   a  normal  population  is  not 
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valid  and  the  data  points   are  not  inaepenaent.    Tnus  tne 
student's  t-test  is  rejecting  H0  when  it  shoula  oe  acceptea. 

4.2    EXAMINATION  QF  SOLAR  RADIATION  MODEL 

Two  approaches  were  usea  to  determine  optimal  values 
for  the  constants,  r  and  X,  in  Garwooa's  soiar  raaiation  pa- 
rameterization. In  the  first  method,  regression  analysis 
was  used  to  determine  the  values  of  r  ana  2f  that  proviaeo 
the  best  agreement  between  Garwood's  ana  Zaneveid's  parame- 
terizations.  Then  r  ana  ^f  were  cnangea  ana  m-  ana  p  were 
varied  to  minimize  E  for  1977. 

Since  Garwooa's  parameterization  is  an  approximacion  to 
a  two  banawidth  moael  of  solar  raaiation  aDsorption,  a  non- 
linear regression  was  performed  to  compare  Zaneveia's  param- 
eterization with  that  of  Paulson  ana  Simpson  (see  Taole  1). 
Zaneveia's  model  of  solar  raaiation  absorption  is  an  empiri- 
cal representation  of  the  extensive  aosorptance  ana  trans- 
raittance  data  set  collectea  by  Spinraa  et.  ai.  (1979);  tnus 
it  is  assumed  that  agreement  with  Zaneveia's  moael  is  tanta- 
mount to  agreement  with  the  data  set.  The  values  wnicn  gave 
the  best  agreement  in  the  upper  10m  were  (r  =0.4,  &J  = 
0.08m")  ana  (r  =  0.6  and  ^  =  3.0m").  The  assumption  tnac 
one  of  the  two  extinction  deptns  is  much  smaller  tnan  tne 
other  is  verified  in  that  #a  is  approximately  two  oraers  of 
magnitude  larger  than  ^  .  The  upper  10m  was  cnosen  since 
the  absorption  of  solar  raaiation  above  the  extinction 
depth,  is  most  critical  to  the  mixea  layer  neat  ana  ouoyancy 
budgets.  Linear  regression  analysis  was  also  usea  to  com- 
pare the  logarithm  of  Garwood's  parameterization  to  tne  log- 
arithm of  Zaneveia's  empirical  formulation  roc  -10m  <  z  < 
-lm.  That  is,  the  equations 
ln(r)  +  X   z 

and 

y  -I 

0  z  +  ln(l  -  a  tan  (bx) ) 

were  compared  to  obtain  the  values  of  r  ana  o  wnicn  gave  the 
best  agreement.  The  best  agreement  was  obtained  for  (r  = 
0.4,  0  =  0.08m"),  wnicn  again  verifies  that  Garwooa's  param- 
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eterization  approximates  the  aoubie  exponential  mooei,  wnich 
was  founo  to  aaequately  descrioe  tne  aosorption  of  soiar  ra- 
diation witn  respect  to  mixed  layer  dynamics  (Simpson  ana 
Dickey,  1981)  . 

We  are  concerned  with  the  effects  of  the  parameteriza- 
tion of  solar  radiation  on  the  neat  ouuget  of  the  mixed  .Lay- 
er. With  this  in  mind,  Zaneveid's  and  Garwood's  parameten- 
zations  were  compared  oy  calculating  the  cnange  in 
temperature  of  an  initially  nomogeneous  water  column  aue  to 
the  aosorption  of  soiar  radiation.  The  initial  temperature 
of  the  water  column  was  12C  ano  the  equation 

a  Z      £Cp    Bz 

was  integrated  for   12  hours.    The  surface   soiar  raaiation 
had  a  cosine   dependence  in  time  v/itn  a  peaK   at  local  noon. 
The  temperature  profiles  after  12  hours  for  three  cnoices  of 
r  ana  o  are  shown  in  Fig.  2a.   The  results  for  (r  =  0.4,  o  = 
0.0  8m  )   and  (r  =  0.45,   0  =  0.0  95m  )   ooth  overestimate  tne 
temperature  at  lm,   underestimate  the   temperature  at  2m  and 
slightly  overestimate  tne   temperature  from  4  to  10m.    Tne 
set  (r  =  0.5,   #=   0.11m  )   estimates  tne  temperature  at  lm 
ano  2m  well  out  overestimates  the  temperature  from  2  to  10m. 
The  change  in  heat  content  is  Dest  estimated  Dy  ( r  =  0.45,  # 
=  0.095m  ),   since   tne  underestimate  of  the   temperature  in 
the  first  2m   is  approximately  baianceu  oy   tne  overestimate 
of  the  2  to  10m  temperatures.     Next,   the  value  of  r,   tne 
fraction  of  solar  radiation  wnicn  penetrates  1  m,   was  neic 
constant  ano  o  was  variec.   For  r  =  0.5  ano  0.45,  0  =  0.11m 
and  o  =  0.095m",   respectively,   gave  tne  Dest  agreement  ue- 
tween  Zaneveid's  and  Garwood's   parameterizations.    In  Fig. 
20,    profiles  are   plotted  for   r  =   0.50  ano  X  =   0.10m  , 
0.15m",  and  0.20m"'.   Fig.  2c  is  similar  to  Fig.   20,  Dut  ror 
r  =  0.45.    In  both  of  these  figures,  all  three  cnoices  or  tf 
give  good  agreement  near  the  surface  wnere  tne  gradient  witn 
depth  is  very  large.    Compared   witn  Zaneveid's  curve,   tne 
curve  for   which  X  =   0.20m   turns  too  soon,    which  implies 
that  the  scale  depth  for  aosorption  is  too  deep.    Tne  curve 
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for  which  0  =  0.10m"  has  approximately  equal  area  aDOve  ana 
below  Zaneveld's  curve;  which  makes  the  integratea  heating 
in  the  upper  10m  the  same. 

The  annual  cycle  for  1977  was  simulated  with  r  =  0.50 
and  %  =  0.10m"  ,  0.15m  and  0.20m  for  a  range  of  values  of  m- 
and  p  .  The  values  of  E  were  higher  for  o*  =  0.15m  ana  V  = 
0.20m*'  than  for  #=  0.10m  and  the  minima  were  displacea  no 
higher  m_  values  for  higher  ^  values.  The  cnange  in  m-  can 
be  explained  since  larger  X  values  imply  snaliower  extinc- 
tion depths.  Therefore,  more  solar  raaiacion  is  aosorbeo 
near  the  surface  and  larger  surface  stresses  are  requirea  to 
mix  the  more  staDle  water  column.  Using  the  bootstrapping 
method  or  the  PPP,  there  is  no  significant  difference  among 
results  for  X  =  0.10m"f  ^=  0.15m"'  and  X  =  0.20m"1  at  the  90% 
confidence  level,  where  (m  =  2.5,  p  =4.0). 

4.3   COMPARISON  WITH  PREVIOUSLY  CALCULATEP  VALUES 

The  optimal  values  of  the  moael  constants  found  in  tms 
study  are  quite  different  from  the  values  previously  used 
and  emphasize  the  need  for  careful  evaluation  of  tne  con- 
stants in  different  conditions.  The  values  of  these  con- 
stants determined  oy  Garwooa  (1976)  were  m-  =  7.6,  p_  =  m^/iUj 
=  1,  r  =0.5  and  O  =  0.1m  .  Tne  Garwooa  model  nas  previ- 
ously been   tested  with  OS  Papa   data  Dy  Garwood   ana  Aaamec 

(1982)  ,  who  specified  m-  =  0.75,  p  =  1.0,  r  =  0.5,   ana  ^  = 
-i  *         3 

0.2m  .  The  ocean  thermal  structure  ana  surface  flux  infor- 
mation used  in  the  present  stuay  is  quite  uifferent  from 
chat  used  in  tne  OS  Papa  simulations.  The  ocean  thermal 
structure  usea  in  the  OS  Papa  cases  was  aerivea  from  XBT  oo- 
servations  taken  at  the  same  site  approximately  every  two 
days.  The  thermal  structure  for  tnis  stuay  was  oDtainea 
from  monthly  objective  analyses  of  XBT  ooservations  on  hori- 
zontal surfaces  of  constant  depth.  Also,  tne  XBT  ooserva- 
tions can  resolve  temperature  graaients  with  ^z  <  5m, 
whereas  the  TRANSPAC  analyses  ao  not  resolve  graaients  witn 
£±z  <  20m.  In  the  OS  Papa  simulations  of  Garwooa  ana  Aaamec 
(1982)  ,   the  surface  stress  and   the  heat  fluxes  were  calcu- 
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lated  from  3-h  ooservations  of  SST,  air  Lemper a cure,  aew 
point  temperature,  cioua  cover,  wina  speea  ana  direction. 
The  v/ma  speed  ana  direction  were  aajuscea  to  10  m.  For 
this  study,  the  surface  stress  is  caicuiatea  from  a  moael 
producea  surface  wind  speea  aajuscea  co  19  m.  The  cocai 
heat  flux  ana  solar  raaiation  flux  are  caicuiacea  from  moael 
heat  ana  moisture  fieias. 

Model  sensicivicy  co  Variations  of  the  vaiues  or  tne 
moael  constants  may  be  juaged  oy  comparing  time  series  of 
SST  generacea  using  different  values  of  che  constants  wich  a 
time  series  of  tne  analyzea  SSTs.  Tne  improvement  in  model 
performance  using  che  values  aeterminea  in  this  stuuy  com- 
pared with  eicher  the  Garwooa  values  of  cne  conscancs  or  cue 
OS  Papa  values  is  illuscracea  in  Figs.  3a,  o  ana  c.  In  Fig. 
3a,  using  values  which  are  cypical  ror  this  study  (m_  =  2.0, 


3 


P^  =  2.0,  r  =  0.5,  X  =  0.1  m  ),  cne  moael  accuracely  simu- 
lates the  annual   cycle  of  sea  surface   temperature.    Using 

the  OS   Papa  values  (m_=   0.75,   P,  =   1.0,   r  -   0.5,    X  = 
-I  3  3 

0.2m  ) ,   the  maximum   summer  temperature  is  coo   hign  uue  co 

too  little  wind  stirring  (m-  =  0.75)  ana  aue  co  cne  snaiiow- 
er  extinction  aepth  (^=  0.2m")  usea  in  aosorpcion  of  cne 
solar  raaiation  (Fig.  3b).  Furcher  there  is  less  variabili- 
ty in  the  summertime  mixea  layer  deptn  ana  cne  race  of  cne 
fall  mixed  layer  deepening  is  less  tnan  in  Fig.  3a.  In  Fig. 
3c,  using  the  values  aeterminea  oy  Garwooa  (m_  =7.6,  p_  = 
1.0,  r  =  0.5,  0  =  0.1m  ) ,  che  moael  aoes  noc  warm  enougn  in 
the  summer.  This  is  aue  to  cne  race  chac  m3  is  coo  large, 
causing  the  effecc  of  the  wina  scirring  co  oe  coo  greac  ana 
keeping  tne  mixed  layer  mixed  too  deep.  Tne  variability  in 
the  summertime  mixed  layer  aepen  ana  cne  race  of  cne  rail 
deepening  is  much  greater  tnan  in  Fig.  3a.  Also  tne  winter 
mixed  layer  is  too  deep  compared  to  Fig.  3a.  Anocner  impor- 
tant difference  among  cne  three  simulations  is  aepen  of  cne 
summer  mixed  layer  ana  the  time  of  cne  occurence  of  cne 
spring  transicion.  The  spring  transicion,  wnicn  is  cne 
abrupt  ana  severe  recreac  of  cne  mixea  layer,  occurs  near 
cne  ena  of  Feoruary  ror   tne  simulacions  using  cne  conscancs 
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determined  in  this  study  and  the  OS  Papa  constants.  How- 
ever, in  the  simulation  using  the  Garwooa  constants,  the 
transition  is  delayed  until  mid-April  ana  the  upper  ocean 
never  achieves  the  shallow  mixed  layer  ana  strongly  strati- 
fied seasonal  thermocline  that  are  cnaracteristic  or  the 
summer  mixed  layer  in  mid-latitudes. 

5.   SUMMARY  AND  CONCLUSIONS 

The  optimal  values  of  the  constants  for  any  turoulence 
model  are  expectea  to  vary  for  different  specifications  or 
the  atmospheric  forcing  ana  the  ocean  thermal  structure.  We 
have  shown  that  the  values  of  the  moael  constants  estimated 
from  laboratory  data,  and  Ocean  Weatner  Snip  ooservations 
yiela  quite  different  results,  wnen  appliea  to  the  same  sim- 
ulation (Figs.  3a,  3b  and  3c) ,  than  tne  "optimal"  values  ae- 
termined  here.  Whether  it  woula  be  Getter  to  alter  tne  val- 
ues of  the  exchange  coefficients  for  the  surrace  riuxes 
rather  than  the  values  of  the  closure  constants  is  largely  a 
philosophical  question.  Knowieage  of  surrace  fluxes  ana 
geophysical  turbulence  is  not  yet  aaequate  to  aistmguisn 
between  the  two  approaches.  We  nave  cnosen  the  seconc  ap- 
proach largely  for  its  inherent  simplicity  ana  to  aemon- 
strate  the  importance  of  understanding  tne  potential  inter- 
actions between  the  model  and  the  aata  usea  to  initialize 
and  force  the  moael. 

The  optimal  values  may  also  cnange  for  different  atmos- 
pheric and  oceanic  regimes,  oecause  ail  possioie  physical 
processes  are  not  parameterizea  in  the  moael.  We  nave  snown 
that  the  Garwooa  moael  is  relatively  insensitive  to  varia- 
tions in  the  model  constants  within  a  fairly  oroaa  region  of 
parameter  space.  However,  the  year-to-year  variations  in 
the  optimal  values  of  the  moael  constants  are  significant. 
In  particular,  1976  was  founa  to  be  significantly  different 
than  1977  ana  1978  by  all  three  statistical  tests  employea 
in  this  study.  This  reflects  the  fact  that  tne  atmospheric 
forcing  and  the  response  of  the  ocean  thermal  structure  were 
quite  different  among  the  three  years  (White,  et.  ai. ,  1980; 
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Haney,  1980;  ana  Elsberry,  1983).  An  improvement  in  moael 
skill  could  be  achieved  by  calculating  the  optimal  values  or 
the  model  constants  for  each  year.  Because  tne  atmospheric 
forcing  and  the  ocean  thermal  structure  are  caicuxatcu  in 
the  same  manner  for  all  three  years,  tne  variations  in  cne 
optimal  values  of  the  constants  prooaoiy  represent  cne  er- 
fects  of  physics  not  incorporated  in  the  moael  ana  noise  in 
the  analyses  of  ocean  thermal  structure  ana  tne  surface 
fluxes.  For  these  reasons,  the  same  vaiues  will  ue  usea  for 
all  experiments  using  this  type  or  forcing  in  the  Worth  Pa- 
cific.  The  values,  p   =  4.0,  m,  =  2.0,  r  =  0.5,  ana  ^  =  0.1 

-I  3         3 

m  ,   were   cnosen  as  oeing   representative  or  tne   swarms  or 

points  from  the  three  years  tnat  ao  not  airrer  signir icantly 
at  the  90%  confiaence  level. 

Three  statistical  tools,  tne  pool  ana  permutation  pro- 
cedure, the  bootstrapping  technique  ana  tne  stuaent's  t- 
test,  were  used  to  aetermine  wnetner  cifferent  parameter 
values  gave  statistically  aifferent  moael  results.  In  gen- 
eral, the  bootstrapping  tecnnique  generatea  siigntiy  oroaaer 
reference  distributions;  however,  it  influencea  tne  infer- 
ences arawn  from  tne  tests  in  only  one  case  (1976:  m-  =3.5 
vs.  m^  =  2.0)  .  The  stuaent's  c-test,  on  tne  o trier  iiana, 
leaa  to  rejection  of  tne  null  hypotnesis  in  far  too  many 
cases.  This  points  out  tne  nazard  of  assuming  Gaussian 
pdf ' s  for  this  type  of  aata.  We  conciuae  tnat  statistical 
tests  which  put  less  weight  on  the  inaepenoence  of  the  aata 
points  are  more  useful  for  stuaies  or  cms  nature. 

The  aoscrption  of  solar  raaiacion  with  depth  is  not 
consiaerea  one  of  the  funaamentai  aspects  of  geophysical 
turoulence  moaeiing.  However,  we  nave  shown  tnat  its  proper 
specification  is  important  for  simulations  of  tne  annual  cy- 
cle and  crucial  to  simulations  of  tne  summer  mixea  layer. 
The  parameterization  employee  in  the  Garwooa  moael  reproauc- 
es  the  results  of  more  complicatea  formulations  ana,  witn 
optimizea  coefficients,  improves  the  summer  ana  annual  simu- 
lation of  tne  upper  ocean  temperature. 
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month.  The  moaei  constants  aeterminea  in  tnis  stuay 
(m,  =  2.0,  p,  =  4.0,  r  =  0.5,  4  =  0.1m*')  were  usea  in 
the  simulation. 

8.  Fig.  3b  Similar  to  Fie.  4a  except  using  OS  Papa  moaei 
constants  (m,  =  0.75,  p3    =  1.0,  r  =  O.o,  a  =  0.2m  ). 

9.  Fig.  3c  Similar  to  Fig.  4a  except  using  Garwooa 
(1976)  model  constants  (m_  =  7.6,  p_  =  1.0,  r  =  O.b, 
*  =  l.Om"1)  .  *         J 
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Fig.  Id  Similar  to  Fig.  la  except  for  1977. 
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Fig.  lc  Similar  to  Fig.  la  except  for  1978 
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Fig.   2a  Temperature  profiles  in  tne  upper  10m  or  the 
in  for   Zaneveia's  empirical  formuiamon   ana  Gar- 
Da's   formulation  using   tnree   r^d  pairs:    (0.40, 


oceai 


0.080),  (0.45,  0.095)  ana  (0.50t  O.llO).  For  eacn 
value  of  r  a  regression  was  perrormea  against  Zane- 
veia's ..empirical  formulation  to  rina  tne  optimal  val- 
ue of  6. 
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Fig.  2b  Temperature  profiles  in  tne  upper  10m  or  tne 
ocean  for^  Zaneveia's  empirical  formulation  ana  ror 
tnree  (r,o)  pairs  of  Garwooa's  parameterization.  In 
tnese,  r  was  set  to  0.50  ana  if  was  variea  (0.10, 
0.15,  ana  0.20)  . 
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Fig.    2c   Similar    to   Fig.    2o    except    ror    r    =    .45 
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